Spin-orbit induced backflow in neutron matter with auxiliary field diffusion Monte 

Carlo 
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The energy per particle of zero-temperature neutron matter is investigated, with particular em- 
phasis on the role of the L ■ S interaction. An analysis of the importance of explicit spin-orbit 
correlations in the description of the system is carried out by the auxiliary field diffusion Monte 
Carlo method. The improved nodal structure of the guiding function, constructed by explicitly con- 
sidering these correlations, lowers the energy. The proposed spin-backflow orbitals can conveniently 
be used also in Green's Function Monte Carlo calculations of light nuclei. 
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I. INTRODUCTION 

In recent investigations 0, Q] of the ground state and 
the magnetic properties of neutron matter with mod- 
ern nuclear interactions of the Urbana-Argonne type 
0] good agreement was observed between results ob- 
tained with the Auxiliary Field Diffusion Monte Carlo 
(AFDMC) 0, a calculation by Morales, Jr. et al. @, 
performed with the Variational Chain Summation (VCS) 
method @, Q, and Brueckner Hartree Fock estimates 

However, in spite of the overall agreement of the equa- 
tion of state and the spin susceptibility, marked differ- 
ences exist between the AFDMC and the VCS calcula- 
tions, concerning the contribution to the energy due to 
the spin-orbit component of the two-body interaction. 
VCS calculations, performed with the Argonne vig Q 
and the Urbana-IX three-body potential (AU18 Hamil- 
tonian), find large and negative contributions from the 
cluster terms Cls with either spin-orbit correlations in 
the trial function and/or the spin-orbit potential. For 
instance, in correspondence with optimal trial functions, 
Cls amounts to -5.8 MeV at a density p = p n = 0.16 
fm~ 3 and -12.1 MeV at twice the same density. The 
AFDMC calculations of Ref. p) , performed with the sim- 
plified, but still realistic, version of the Argonne two- 
body potential, vs> 3, plus the Urbana IX three-body 
interaction (AU8' Hamiltonian) get energy differences 
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AEr,s between the AU8' Hamiltonians with and without 
the spin-orbit potential which are small and positive. At 
p = p a , AE LS = 0.2 MeV and at p = 2p , AE LS = 1.12 
MeV. 

Such a discrepancy cannot be ascribed to differences 
between the two potentials v\s and vg> . It is well known 
that they provide very close results for the energy per 
particle Q . In addition, this discrepancy is confirmed by 
other FHNC/SOC calculations pj, performed by using 
the approximations of Ref. and the AU8' Hamilto- 
nian. They give AE LS = -3.7 MeV at p and AE LS = 
-10.1 MeV at 2p , with C LS being -3.9 MeV and -10.7 
MeV respectively. 

A possible source for this disagreement might be the 
use of a less than satisfactory guiding function in the 
AFDMC method The nodes of the plane wave Slater 
determinant might be too poor, particularly when the 
interaction includes a spin-orbit potential, like in At/8'. 
The results of the FHNC/SOC calculations of Ref. Q 
with the AU8' Hamiltonian and a trial function of the 
type Fq (not containing spin-orbit correlations), give val- 
ues for AEls quite close to the AFDMC ones, which may 
confirm this hypothesis. 

To clarify this issue we modify the guiding functions in 
our quantum Monte Carlo calculations to contain explicit 
spin-orbit correlations. This is efficiently done by con- 
sidering orbitals of the spin-backflow form in the Slater 
determinant, as explained below. 

The structure of this paper is as follows. In Sec. [H] 
we show the guiding function used. The computational 
details are given in Sec. IIIII The results are shown and 
discussed in Sec. IIVI The conclusions and perspectives 
of the present work can be found in Sec. 



II. SPIN-ORBIT INDUCED BACKFLOW 

The L- S correlation in FHNC/SOC and VCS calcu- 
lations takes the form 
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F b (l, 2) = -f b (r 12 ) [r 12 x (Vi - V 2 )] • (a 1 + a 2 ) (1) 

4i 

where ryi = \ v n\ and &j are the Pauli matrices for the 
jth particle. 

By inspection of the cluster terms of (AU8 1 ) at the 
two-body order, one observes that the leading of the 
spin-orbit terms are those usually denoted as bbc and 
ebb in FHNC/SOC theory 12]. In the calculations of 
Ref. these terms are responsible for ~ 80% of the to- 

(2) 

tal contribution from all the C L g terms in the density 
range (3/4)p < P < (5/2)po- One can easily prove that 
exactly the same expressions of the ebb and bbc terms 
are obtained by the following simplified L ■ S correlation. 



F b (l,2) = -/ 6 (ri 2 )[r 12 



(Vi x o-i - V 2 x a 2 )\ ■ (2) 



It is found that, at the two-body level of the FHNC/SOC 
theory, F), leads to an energy which is only ~ 10% differ- 
ent from that obtained with Ff,. 

The important feature of the L ■ S correlation of Eq. 
(J2J) for AFDMC calculations is that, similarly to the case 
of standard backflow , it can be implemented in quan- 
tum Monte Carlo simulations by substituting the plane 
wave orbitals of the Slater function with the following 
spin-backflow ones 
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FIG. 1: Neutron matter energy per particle in MeV as a func- 
tion of the strength (3 of the spin-orbit correlation function 
at p — 2po for 14 neutrons. Solid line stands for the AU6' 
interaction, dotted line for the AU8' potential. 



where m is assumed to be the average of the neutron 
and proton masses and h 2 /m = 41.47108 MeV fm 2 ; the 
two-body and three-body potentials, Vij and Vjki, are 
the Argonne v' s and the the Urbana IX potentials Q. 
The three-body contributions to the energy are large, 
particularly at high densities, and cannot be neglected 
even in a survey calculation like ours. 

In neutron matter this interaction can be written in 
terms of four components: 



0fcC?') = eX P ik ' r 3 



p x - 
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h{r jk ){r jk x fc) 



where /? is a spin-orbit strength parameter. For [3 = 1 
there is a direct correspondence of with Ft,, This 

spin-backflow ansatz can also be used in the GFMC sim- 
ulations of small nuclear systems 0, ITil ITHj l , to include 
L ■ S correlations. 

We present and discuss, in the following, the results 
obtained in AFDMC simulations of neutron matter en- 
ergy with the AU8' Hamiltonian and the nodal surface 
of the spin-backflow Slater function. We will show that 
these nodes serve to lower the AFDMC energy per parti- 
cle of neutron matter by a sizable amount, which however 
is too small to solve the CVS and AFDMC discrepancy, 
particularly at higher densities. 



III. CALCULATION DETAILS 

We describe the neutron system by the non relativistic 
Hamiltonian 



H = 



T + V2+V3 
2m 



= -iE V 2 +5> 3 ,+ £ V jkl , (4) 



OP(ij) = 1,^-^,^,1,-S 



(5) 
(0) 



where Sij and L ■ S arc the usual tensor and spin-orbit 
operators. The functions v p {rij) can be found on Ref. 
and also on Ref. pj. 

The Urbana-IX three-body interaction is given by the 
sum 



V jk i 



ySI 



ySD 
v jkl 1 



(7) 



j=l,N 



3<k<l 



where Vj k j is a spin independent three-body short range 
part, and the spin-dependent part, V^ff, in neutron mat- 
ter, reduces to a sum of terms containing only two-body 
spin operators, with a form and strength that depends on 
the positions of three particles. Their explicit expressions 
can be found in Ref. [l|. 

We have also considered an interaction obtained from 
AU8' by dropping the spin-orbit term, and, as in Ref. 0, 
it is denoted as AU6'. 

The AFDMC 5] method used in these calculations has 
been previously described in detail 0. It allows Monte 
Carlo simulations to be performed on a relatively large 
nuclear system at the required accuracy, thanks to the 
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introduction of auxiliary fields, which uncouple the spin- 
dependent interaction between particles by means of a 
Hubbard-Stratonovich transformation. While the prop- 
agation of the particle coordinates is done as in diffu- 
sion Monte Carlo, that of the spin variables, after having 
sampled the auxiliary fields, results in a rotation of each 
particle's spinor. 

The guiding function used in this work is given by a 
Jastrow product correlating an TV x N Slater determinant 

(R, S\* JSB ) = [] f(r jk )$ B (R, S), (8) 
j<k 

where R and S denote the set of particle coordinates and 
of spinors respectively. The space and spin orbitals in 
$>b(R,S) are given in Eq. ©, and they operate on the 
spin states in the following way 



ik -Ti 



(cosh(^ fc «)+i|(*)sinh(A fc (*))) | t) 
(i|(i) + ayi))sinh(^ fe (i))||>] , (9) 



where 



(10) 



and {a — x,y,z} are the components of an unit 

vector in the direction of Au(i). A similar expression 
can be obtained for the action over the spin down single 
particle state. We denote this function by JSB as op- 
posed to JS which stands for a function without explicit 
backflow correlations, i. c, with simple plane waves in 
the spatial part of the orbitals in the Slater determinant 
(f b (r) = 0). The AFDMC method described in Ref. Q 
needs only slight modifications to deal with {R, S\^jsb) 
as guiding function. They refer mainly to calculating the 
kinetic energy part of the local energy, and the gradient 
of the guiding function in the drift of the walker. 

The Jastrow and spin-orbit correlation functions f(r) 
and fb(r) have been taken as the first and fourth com- 
ponents respectively of the FHNC/SOC correlation op- 
erator that minimizes the energy per particle of neutron 
matter at the desired density [TO ]. 



IV. RESULTS 

We have made our calculations within the full simula- 
tion box, and we have taken into account the 26 neighbor- 
ing boxes in the tabulation of the correlations f(r) and 
fi(r), and of the various components v p (r) of the two- 
body potential, as described in Ref. pj. Our results are 
therefore already tail corrected for a Hamiltonian with 



TABLE I: AFDMC energies per particle in MeV for the AU6' 
and AU8' interactions obtained for a system of 14 or 66 neu- 
trons in a periodical box with the guiding functions, ^, JS 
and JSB (with (3 = 1), at po- Error bars for the last digit are 
shown in parentheses. 





N = 14 


N = 66 


JS(AU6') 


48.27(9) 


53.11(9) 


JS(AU8') 


48.4(1) 


54.4(6) 


JSB(AU8') 


46.8(1) 


52.9(2) 



two-body force only. Tail corrections for the three-body 
potential were not included. However previous analyses 
[lj have shown that they are small for systems with 66 
neutrons. 

Fig. shows the AFDMC energies per particle of 14 
neutrons in a periodic simulation cell at p = 2/?o for the 
AU& and AU8' interactions, as a function of the strength 
parameter (3. One can see that, for the AU8' interaction 
the energy minimum is around (3=1, consistent with the 
conjecture that fb(r) moves the nodes in the optimal way 
when the spin-orbit potential is included in the Hamil- 
tonian. We have obtained a lowering of the energy of 1.6 
MeV with respect to the case of the JS nodal surface, 
which corresponds to (3 = in the figure. If we switch off 
the spin-orbit potential by considering the AUG' inter- 
action, the minimum is found at (3 = 0, confirming that 
the spin-backflow nodes are energetically advantageous 
only in the presence of the spin-orbit component of vs> ■ 

AFDMC simulations for the AU8' Hamiltonian with 
the JSB guiding function have been also carried out for 
66 neutrons. As in the case of the 14 neutron system, we 
have found a minimum of the energy at (3 ~ 1 . The result 
obtained is compared in Table [I] with those at (3 = for 
both the AU8' and the AU6' interactions. In spite of 
sizable differences between the energies per particle of 
the 14 neutron and the 66 neutron systems, the gain in 
energy (E(JSB) — E(JS))/N is roughly independent on 
the number of particles in the box. The large differences 
between E(14) and E(66) are mainly due to the effect of 
three-body interaction. It has been shown p| that the 
finite size effects on E(66) are rather small. 

The dependence of the JSB energy on the density is 
reported in Table [n] for the 14 neutron system. There 
is a very weak dependence of AE^g on the density, in 
contrast with VCS results. 

In order to make comparison with recent quantum 
Monte Carlo calculations performed for 14 neutrons in- 
teracting via the v' s two-body potential by Carlson et al. 
[TH , we report in Table IIIII the corresponding AFDMC 
results. The table displays the energy difference AEls 
between the energy obtained with v' s and v' e two-body 
potential. 

The calculations of Ref. to which the results re- 
ported in Table ITTll refer, have been performed with the 
vg' potential cut off at the edge of the box. There- 
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TABLE II: AFDMC energies per particle in MeV for the 
AUG' and AU8' interactions obtained for a system of 14 neu- 
trons in a periodical box with the guiding functions, ^, JS 
and JSB (with (3 = 1), as a function of the density p. Error 
bars for the last digit are shown in parentheses. 



p 


AUG' 


JS 


JSB 


PO 


19.73(5) 


19.76(6) 


18.76(5) 


2po 


48.27(9) 


48.4(1) 


46.8(1) 



TABLE III: Spin-orbit contribution to the energy per parti- 
cle in MeV of neutron matter at density po- The constrained 
(CP) and unconstrained (UC) GFMC results, as well as the 
VCS ones are taken from Ref. 0. The AFDMC results ob- 
tained with the JS guiding function are taken from Ref. Q. 
Error bars for the last digit are shown in parentheses. 

method AEls 

GFMC-CP -1.26(4) 

GFMC-UC -2.9(3) 

AFDMC-JS -0.14(6) 

AFDMC-JSB -1-2(1) 
VCS -3.8 



strained path approximation (GFMC-CP). The VCS es- 
timate seems to be too large. 



V. CONCLUSIONS 

In this paper we have proposed a new kind of space- 
spin orbitals with a spin-backflow form, which is partic- 
ularly useful for taking into account the spin-orbit in- 
teraction of nuclear systems in quantum Monte Carlo 
simulations. An efficient parameterization of the spin- 
backflow is obtained from the spin-orbit correlation of 
FHNC/SOC theory. These spin-backflow orbitals can be 
conveniently used also in other quantum Monte Carlo cal- 
culations, for instance the Green's Function Monte Carlo 
simulations of small nucleon systems. The nodal surface 
provided by this new guiding function is able to decrease 
the energy about 5 per cent. This amount is not suf- 
ficient to solve the spin-orbit discrepancy between the 
variational chain summation and the AFDMC results, 
however the AFDMC results are in good agreement with 
constrained GFMC simulations. 
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fore, the values of AE^s extracted from there might not 
be completely comparable with ours, since ours are al- 
ready tail corrected and have been obtained without in- 
troducing any discontinuity in the potential. AFDMC 
seems to agree reasonably well with GFMC in the con- 
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